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Abstract 

We prove the existence of self-similar solutions to the Fradkov model for two- 

^~~^ dimensional grain growth, which consists of an infinite number of nonlocally coupled 

' ' transport equations for the number densities of grains with given area and number of 

Z^ neighbours (topological class). For the proof we introduce a finite maximal topological 

, class and study an appropriate upwind-discretization of the time dependent problem in 

^ self-similar variables. We first show that the resulting finite dimensional differential sys- 

^"5 tem has nontrivial steady states. Afterwards we let the discretization parameter tend 

l/^ to zero and prove that the steady states converge to a compactly supported self-similar 

^"^ solution for a Fradkov model with finitely many equations. In a third step we let the 

^_^ maximal topology class tend to infinity and obtain self-similar solutions to the original 

p I system that decay exponentially. Finally, we use the upwind discretization to compute 

.^ self-similar solutions numerically. 

^' 
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1 Introduction 

Grain growth denotes the late stage coarsening of polycrystahine materials when migration 
of grain boundaries due to capillary forces causes small grains to vanish and larger grains to 
grow. One often observes that despite having different histories many materials eventually 
exhibit universal statistically self-similar coarsening behaviour, usually referred to as normal 
grain growth. Different approaches have been used to predict and explain this phenomenon, 
see [2l[3] for Monte-Carlo methods, [H] for a study of vertex models, and more recently (TJE] 
for boundary tracking methods. We also refer to the theoretical approach in [5; "5] , which is 
based on the grain boundary character distribution. 

However, it remains a challenge to establish such universal long-time asymptotics in math- 
ematical models and it is often difficult to prove only the existence of self-similar solutions. In 
this article we investigate the existence of self-similar solutions to a kinetic mean-field model 
that has been suggested by Fradkov [lOj to describe grain growth in two dimensions. 

1.1 Pradkov's mean-field model 

We briefly describe the derivation of Fradkov's model and refer to [10 1 111 ^ 112] for more details. 
Our starting point are two-dimensional periodic networks of grain boundaries that meet in 
triple junctions (Fig. [I|. In the case of constant surface energy and infinite mobility of triple 
junctions, the grain boundaries move according to the mean curvature flow while all angles at 
the triple junctions are 27r/3. In this setting one can easily derive the von Neumann-Mullins 




Figure 1: Cartoon of a 2D network of grains with triple junctions, where the Herring condition 
implies that all angles are equal to 27r/3. The encircled numbers refer to the topology classes. 

law for the area a{t) at time t > of a single grain with n edges |17j : 

^a(t) = Ma^(n-6). (1) 

Here M denotes the mobility of the grain boundaries and a the surface tension. 




Figure 2: Neighbour switching 

The evolution of such a network by mean curvature is well-defined \15\ [T6] until two 
vertices on a grain boundary collide, after which topological rearrangements may take place. 



If an edge vanishes an unstable fourfold vertex is produced, which immediately splits up again 
such that two new vertices are connected by a new edge. As a consequence two neighbouring 
grains decrease their topological class (i.e., the number of edges), whereas the other two 
grains increase it (Fig. ^. Furthermore, grains can vanish such that some vertices and 
edges disappear. Due to the von Neumann-Mullins law this can only happen for grains with 
topological class 2 < n < 5. As illustrated in Fig. [3j the vanishing of a grain of topological 
class n = 4 or n = 5 can result in topologically different configurations. 







Figure 3: Grain vanishing 

In order to derive a kinetic description of the evolution of the grain boundary network we 
introduce the number densities /„ (a, t) of grains with topological class n > 2 and area a > 
at time t >0. As long as no topological rearrangements take place, the von Neumann-Mullins 
law (IT]) implies that /„ evolves according to 

dtfn {a, t) + {n- 6) dafn {a, t) = . 

This equation needs to be supplemented with boundary conditions at a = for n > 6. It 
is reasonable to assume that no new grains are created during the coarsening process, which 
implies that 

/„ (0, t) = for n > 7 . (2) 

To model the topological changes we define a collision operator J that couples the equations 
for different topological classes. More precisely, we introduce topological fluxes r/+ and rj~ 
that describe the flux from class n to n + 1 and from n to n — 1, respectively, and set 



iJf)n = r?+_i + Vn 



+1 



Vn 



Vv 



with r]f = 7]2 =0 due to n > 2. Employing a mean- field assumption Fradkov (TU] suggests 
that the fiuxes are given by 



Vn = r/3 nfn 



Vn = r(/3 + l)n/„, 



(3) 



where the coupling weight F describes the intensity of topological changes and depends on the 
complete state of the system in a self-consistent way, see ([6]) below. Moreover, the parameter 
/3 measures the ratio between switching events and vanishing events. Our analysis requires 
/? G (0, 2), but the numerical simulations work well also for larger /3. 

Although Fradkov's model has no upper bound for the topology class n, it is convenient 
for the mathematical analysis to consider variants of the model with 2 < n < N for some 
6 < A'' < oo. In this case we close the equations for {fn)2<n<N by assuming 'r]N_^_i = rj^ = 0. 

Assumption ([3]) implies that the collision terms are given by Jf = TJf with 

( J/)2 = 3 (/3 + 1) /3 - 2/3/2, 

(^/)n = (/3 + 1) (^ + 1) /n+i - (2/3 + l)n/„ + /3(n -!)/„_! for 2 < n < A^ - 1, (4) 

(J/)^ = (3{N-l)fN_i - (/3 + l)NfN, 



where the last identify is not used for N = oo. Notice that this definition ensures the zero 
balance property 



TV 



Y^ (Jf)^ (a, t) = for ah a, t > 0, 



n=2 

which reflects that the number of grains with given area does not change due to switching or 
vanishing events. To summarize, the kinetic model we consider in this paper is given by 

dtfn (a, t) + {n- 6) daU («, t) = T{f{t)) {Jf)^ (a, t) , (a, t) G (0, oo)^ , n > 2, (5) 

with boundary conditions pi), either N < oo or N = c«, and {Jf)n given by Q. 

It remains to determine the coupling weight T in dependence of /. The key idea is to choose 
r such that the total area 



A(i) = Vyn(i) with Yn{t)= afn{a,t) 



is conserved during the evolution. One easily checks that dA/dt = P, where P is the polyhe- 
dral defect defined by 



1^ /-oo 

P{t) = 5" (^ - 6)^n(t) with Xn{t) = / fn{a, t) 



The polyhedral formula P = resembles Euler's formula for networks with triple junctions 
and states that the average number of neighbours per grain is 6. We now readily verify that 
dP/dt = holds if and only if 

E(n-6)V„(0,t) 



n=2 

/3NXN{t) - 2 (/? + 1) X2{t) + J2 nXn{t) 



r(/(t)) = "^ ^ , (6) 



n=2 



where we use the convention NXj^ = for A^ = oo. In particular, ([6| guarantees the 
conservation of area and the polyhedral formula provided that the initial data satisfy P = 0. 
We finally mention that well-posedness of Fradkov's model, both for N < oo and N = oo, 
has been established in [12] for /3 G (0, 2). A similar model with finite A^ has been considered 
in [6]. 

1.2 Self-similar solutions and main result 

Self-similar solutions to ([5| take the form 

fn{a,t) = ^^ , S. = ->0, 

where the sequence g = {gn)n>2 of self-similar profiles satisfies 

-2gr,-[i + Q-n)g'^ = T{Jg)^ (7) 



for some positive constant T as well as the boundary conditions 5n(0) = for n > 6. With 
some abuse of notation, we define the moments 

Xn= QniOdi and K = / ^gn{0 d^ 

Jo Jo 

N N 

P = ^ (n - 6)X„ and ^ = Y1 ^n 



and refer to 



as the polyhedral defect and the area of a self-similar solution, respectively. 

It is important to note that each sufficiently integrable solution to M satisfies the analogue 
of the polyhedral formula, and that the coupling weight F depends on the QnS in a self- 
consistent manner. In fact, multiplying ([7]) by ^, integrating with respect to ^, and summing 
over n, we find that the zero balance property of J implies -P = 0. Similarly, if we multiply by 1 
instead of .^, we easily derive the analogue to Q, that means we have T{g) = rnum(5)/rdcn(fi') 
with 

5 N 

rnum(5) = J^(n-6)'5n(0), TC^) = /SA^Xjv " 2 (/3 + 1) ^2 + ^ nX„ . (8) 

n=2 n=2 

The main mathematical difficulty in the existence proof for self-similar solutions stems from 
the fact that the ordinary differential equation ([T]) is singular at ^ = n — 6 and has different 
transport directions for ^ < n — 6 and ^ > re — 6. In this paper we prove the existence of 
weak self-similar solutions for both N < oo and N = oo, where weak solution means that 
each function gn satisfies 

/ g„((e + 6-n)</>'-0)de + (6-re)+5n(OMO)=r(5) / {Jg)jd^ (9) 

Jo Jo 

for all smooth test functions (/) with compact support in [0, cxd). 

Our existence result can be summarized as follows. 

Theorem 1. Let /? S (0, 2) and assume that either 6 < N < oo or N = oo. Then, there exists 
a weak self-similar solution to the Fradkov model that is nontrivial and nonnegative with finite 
area, and satisfies 



N 

E( 



POO 

e^"X„ + y e^^gn{Odq<oo 



for a// < A < In (1 + 1//3). Moreover, for N < oo all functions gn are supported in [0, A^ — 6]. 
These main assertions can be supplemented by the following remarks. 

1. Since T(g) depends on g homogeneously of order 0, the set of self-similar solutions is 
invariant under scalings gn "^ \gn with A > 0. Therefore we can normalize self-similar 
solutions by prescribing the area. 

2. The weak formulation combined with the integrability condition gn € L^(0,cxd) implies 
regularity results. Specifically, in a first step we find that each function gn is continuous 
at all points $, G [0, oo) \ {re — 6}. Using this we then easily show in a second step that 
gn is even continuously differentiable at all points ^ G [0, oo) \ {re — 7, re — 6, re — 5}. 
We also establish further regularity results that characterize the behaviour of gn near 



^ = re — 6, see Lemma 19, and find for large re that gn is continuous also at ^ = re — 6. 



3. The moment estimates from Theorem [T] imply that Qn decays exponentiaUy in ^. More 
precisely, multiplying ([7|) with e'''^ and integrating over [^, oo) gives gn{^) < Cne^^ for 
some constant Cn and all S,- 

4. Numerical simulations as described in Section |4] indicate, at least for N < oo, that for 
each /3 there exists a unique self-similar solution with prescribed area, but we are not 
able to prove this. 

Our strategy for proving Theorem [T] is inspired by the existence proof for self-similar solutions 
to coagulation equations in [9]- In Sectionl2]we introduce a finite-dimensional dynamical model 
that can be regarded as a semi-discrete upwind scheme for (pi) in self-similar variables, and 
involves the discretization length < e <C 1. To derive this scheme we assume that N < oo, 
restrict the rescaled area variable ^ to a finite domain [0, L] with sufficiently large L, and 
impose artificial Dirichlet conditions at .^ = L. Moreover, we identify a discrete analogue to 
([6]) that ensures the conservation of both area and polyhedral formula. Standard results from 
the theory of dynamical systems then imply the existence of nontrivial steady states for each 
sufficiently small e. 

Afterwards we show that these steady states converge as e —t- to a self-similar profile for 
the Fradkov model for N < oo. The proof of this assertion combines two different arguments: 



First, in Section 3.1 we derive suitable a priori estimates that allow us to extract convergent 



subsequences whose limits provide candidates for the self-similar profiles. Second, in Section 



3.2 we analyse the behaviour near the singular points ^ = n — 6 in order to rule out that Dirac 



masses appear in the limit. 



In Section 3.3 we establish the exponential decay of Xn and derive uniform estimates for 
higher moments. The resulting tightness estimates then enable us to pass to the limit A^ — )• oo 
in Section 3.4 Finally, in Section [4] we illustrate that the upwind discretization of ([5]) in self- 



similar variables, combined with explicit Euler steps for the time discretization, provides a 
convenient algorithm for the numerical computations of self-similar solutions. 

2 The discrete dynamical model 

In order to prove the existence of self-similar solutions we study an upwind finite-difference 
discretization of the time-dependent problem in self-similar variables. To that aim we restrict 
the rescaled area variable ^ to a finite interval [0, L] with L G N and L > N — 6, and for each 
K £ NL we consider the grid points ^k = kL/K, such that 

4+1 - 4 = L/K =: e. 

Notice that for each n > 6 the critical area ^ = n — 6 corresponds to one of the grid points, 
that means for each n > 2 we have 

Ckn=n- 6, kn := (n - 6)K/L 

with kg = 0, kn < for n = 2. ..5, and < kn < K for n = 7...N. 
Using the difference operators V^ and V+ with 






e e 

we mimic the transport term — (^ + 6 — n)g'n{C) by the upwind discretization 

-fe + 6 - n)+V+5^ + (6 + 6 - n)_V-g^n + ^tol (10) 

6 



Here 5^" is the usual Kronecker delta and x± denotes the positive and negative part of x, 
that means x± = maxjitx, 0} > and x = x_|_ — x_. 



At a first glance, the Kronecker delta in ( 10 ) seems to be quite artificial, but it is naturally 



related to the singularity of the transport operator. More precisely, for a continuous variable 
^ one easily shows that 

holds in the sense of distributions, where 5n-Q{C) is the Dirac distribution supported in ^ = 
n — 6. Our discretization of the transport operator satisfies a similar identity which can be 
seen by setting g^ = e~^5^ in (10). The Kronecker delta in (10) therefore guarantees that 
the resulting discrete scheme resembles the continuous dynamics even if mass is concentrated 
near the singularities. 



With ( 10 ) the discrete dynamical model reads 



dt 



9n 



- „k 



25n - (a + 6 - n)+V+5n + (6 + 6 - n)_V-g':, + 5f 5n = nJoln 



(11) 



where the coupling weight T will be defined in Section [2. 1[ To close the system ( |11[ ) we impose 
the boundary conditions 



a'n 







for 7 < n < iV, 



^ 







for 2 < n < A^, 



fl2) 



so (11) becomes an evolution equation for the variables g^ with n = 2...A^ and k = 1...K — 1. 
Notice that the boundary conditions for k = stem naturally from (pi), whereas those for 
k = K reflect the cut off in ^. 

2.1 Moment balances and choice of the coupHng weight 

hi complete analogy to the discussion in Section [T] we choose the discrete coupling coefficient 
r such that (11) with (12) conserves the area. In order to identify the correct formula we 



start with an auxiliary result for a discrete moment Z„ with 

K-l 

^n '■= £ 2_^ l^ndm 
k=l 

where /i„ are arbitrary moment coefficients. 
Lemma 2. We have 

d ^"^ 

— Zn - 2Z„ + Cn = TeY, /^n(^/)„, 

where 

K-l 

Cn ■■= (6 - n)+ fi^gl + e XI '^n9n, 



k=l 



(13) 



k=l 





( V+[^^(6 + 6- 


-n)] 


for k < kn, 


Vt:= 


< ^n" 




for k = kn, 




I V-[^^(6 + 6- 


-n)] 


for k > kn- 



(14) 



Proof. Multiplying (11) by e//^, and summing over k = 1...K — 1, give (13) with Cn = Cn + 
Cn + e/^n"9n" and 



K-l 



K-l 



c+ 

Sn 



-£ Y. ^niik + 6 - n)^V+gt Cn ■=^Y1 ^'niCk + 6 - n)_V~gt 



k=l 



k=l 



We now reformulate C„ and C,^ by means of the discrete integration by parts formula 
K2 K2 

k=Ki k=Ki 

To this end we consider the following two cases: 

Case I: 2 < n < 6. By definition, we have ^fc + 6 — n > for all k = 1...K — 1. This 
implies C,^ = and hence 

K-l K-l 

C = -e 5] Mn(6 + 6 - n)V+<7' = (6 - n)fi^^gi + e ^ 5nV-[/x^(6 + 6 - n)], 
fc=i fc=i 

where we used ^0 = ^-nd the boundary condition g^ = 0. 

Case II: 7 < n < N. Here we have 1 < k^ < K — 1 with £^k„ + 6 — n = 0, and hence 

K-l K~l 

C = -e E ;"n(6 + 6-n)V+<7^ = e ^ 5^v-[^^,(efe + 6 - n)]. 

Similarly, we find 

Cn=-^Y1 ^'niik + 6 - n)V-gi = e ^ <?^V+[/x^(4 + 6 - n)] 
fc=i fc=i 

thanks to the boundary conditions g^ = g^ = 0. D 

We next summarize some elementary properties of the coupling operator J . 

Lemma 3. The coupling matrix J satisfies 

N N-1 N 

Yl ^n{Jf)n = Y. (^-+1 - ^n)Pnfn " J^ (^n " ^n-l)(/3 + l)n/„, 
n=2 n=2 n=3 

where 9n denote arbitrary weights. In particular, we have 

N N N 

E (-^Z)- = 0' E (6 - n){Jf)^ = miN - /32/2 + Y^fn. (15) 

n=2 n=2 n=3 

Proof. The definitions from Q imply 

Af JV-1 N 

E ^"("^Z)" = E ^"(/^ + !)(" + l)/n+l + E ^"^(^ - !)/"-! 
n=2 n=2 n=3 

TV-l AT 



E^"/3n/n-E^«(/5 + !)"/«' 



n=2 n=3 

and all claims follow immediately by direct computations. D 

For the following considerations we introduce the discrete moments 

K-l K-l 

X„ := e ^ ffn, 1; := e E ^'^^^ 

fe=i fc=i 



as well as the auxiliary quantity 



N K-1 



Q ■=^^Y1 sgn(Cfc + 6 



n)9n- 



n=2 fc=l 

We also define the discrete area and the discrete polyhedral defect by 

N N 



A:=Y,yn, P:=^(n-6)X„. 



n=2 



Tl=2 



Corollary 4. We have 
d 



dt 

dt 



Xn = Xn+T{JX)^-{6-n)+gl, 



K-l 



Yn = {n- 6)Xn + r( jy)„ + e^Yl ^S^^^^ + ^ 



n)9n 



(16) 
(17) 



fc=i 



Proof. The claim for X„ follows from Lemma ^ since ^„ = 1 implies 



^n 



V+[Cfc + 6-n] = 1 for k < K, 
1 for k = kn, 

V-[Ck + Q-n] = 1 for k > kn- 



Similarly, with /x^ = ^^ we find 

V+MCk + 6-n)] = 2^k + 6-n + e for k < kn, 
& = 2,ffc + 6 - n for A; = kn, 

^~[Ck{Ck + Q-n)] = 2^fc + 6-n-e for k > kn- 

This means 7/„ = 2^fe + 6 — n — esgn(^fc + 6 — n) for all n = 2...N and k = 1...K — 1, so Lemma 
[I implies ([17|. D 



We now show that the initial value problem for the differential system (11) has a global 
unique solution with state space 

U = [{g^^ : A = l, P + eQ = o], 

{\ '^Jn=2...N,k=l...K-l ^ j ^ 

provided that F = rnum/rden is given by 

5 

Tnum := ^ (6 - n)(6 - n + e)^.^, 

n=2 
N N K-l 



(18) 



n=2 



n=2 fc=l 



Lemma 5. Definition ( |18[ ) implies thatU is invariant under the flow of (11) wii/i (12). 
Proof. From Corollary H and (15)i we infer that 

d d ^ ^ 

_^ = P + eQ, _p_p = ^(6-n)25i-F^(6-n)(JX)„. 

n=2 n=2 



(19) 



9 



We now compute dQ/ dt using Lemma [2j With ^„ = sgn(^fc + 6 — n) we find 

-V+[CA: + 6-n] = -1 for A; < /c„, 
Vn= ^ for A; = /c„, 

+V" [^fc + 6 - n] = 1 for A; > A;„ 



and hence /x^ = ry,„. Using (13) and (14) we therefore conclude that 
, 5 N K-l 



dt 



Q-Q = -Y^{Q- n)gl + eF J^ ^ sgn(a + 6 - n) {jg^)^. 



n=2 



n=2 k=l 



Combining this with (19) we get 



d2 d 

-^A =—{P + eQ)-{P + eQ) = Tnum - T • Tden, 



where Tnum and Fden are defined in (18). In particular, V = T^^^/T^^^ imphes d{P+£Q)/ dt = 
dA/ dt = for aU states from lA. D 



2.2 Global solutions and steady states 

Our next goal is to establish the existence of global in time solutions to the discrete dynamical 



system (11). In particular, we show that the restriction /3 < 2 implies that the denominator 



of r is strictly positive for all times. 

Lemma 6. There exist constants Dj^ and en < 1 that depend only on N such that for all 
states in 



U. 



{(5^)GW:5n>0 for n = 2...N, k = 1...K - l] , 



we have 



N 



N 



{2-P)Y,Xn< Tden < I^iV J^ ^. 



(20) 



n=2 



n=2 



and 



2-/3 ,^ ^Dn 

r S i den S : 

L e 



r < 

i num ^ 



16 + 4e 



e^ 



(21) 



provided that < e < e^y. 



Proof. For all states from lA^ we find, using (15)2 and (18), that 

N N K-l 



Lden = -2(/3 + 1)^2 + mXN + 6 ^ ^n " e' ^ ^ sgn(a + 6 - n) (5'= + j/)„, (22) 



n=2 



n=2 fe=l 



where we used that Yln=2 '"'■^n = 6 ^n=2 ^" ~ ^Q thanks to P + eQ = Q. With ( 22 ) and 



K-l 



£2 ^ sgn(efc + 6 - n) (/ + j/), 



fc=i 



K-l 



< e'^Dn ^9^ = eD^Xn 



k=l 



10 



we find 



TV 



N 



(4 - 2/? - eDn) Y,^n< Tden < (/3iV + 6 + eDn) Y. X„, 



n=2 



n=2 



which then implies (20) for all sufficiently small e. In view oi e < ^^ < L — e for all k 



1...K — 1 we have that eX^ <Yn< LX„, 2 < n < N, and conclude that 

N N 

eYXn<l<LY,Xn, 



n=2 



n=2 



SO pO| ) implies (21 )i. Moreover, (21)2 holds because we have Fnum < (16 + 4e) ^^^2 5n ^^'^ 
EL29i < e-'. " D 

Now we are able to prove that the initial value problem for the discrete model is globally 
well-posed with state space W+ . 

Lemma 7. Let Q < e < en- Then, for any initial data from U^ there exists a unique global 



solution to (11)-(12) that takes values in U+ for all times t >0. 



Proof. The estimates from Lemma |6] guarantee that the mapping g £ U i-^ T is locally 
Lipschitz, so local existence and uniqueness of a solution with values in Vl follow from standard 
results. Moreover, due to the upwind discretization of the transport operator we easily show 
that the flow preserves the nonnegativity of 5. Finally, conservation of area and the estimates 
from Lemma [6] imply that T is uniformly bounded in time, and hence the global existence of 
solutions. n 

Since the set U^ is convex and compact, the existence of steady state solutions follows 
from standard results. 



Corollary 8. For all sufficiently small e there exists a steady state solution g G U^ to (11)- 



(12). 



Proof. See, for instance, Proposition 22.13 in [T]. 



D 



Remark 9. If g a steady state solution to (11) and (12) then so is Xg for any A > with the 
same coupling weight T. 

We conclude with further properties of steady state solutions. 



Lemma 10. Each steady state solution g to (11) and (12) satisfies 

5 N 



Y,{Q-n)gl = Y,X^, 



(23) 



n=2 



n=2 



and thus we have T < 



2-/3 



for all E < E]y. 



Proof. Equation p,3v follows from summing over n = 2...A^ in the stationary version of (16) 



with the help of (|T5|). Using (18) and en < 1 we then derive 

5 N 



Fnum < 5 ^ (n - 6)gl, = 5 ^ X„, 



n=2 



n=2 



and combining this with (20) we find the desired result. 



D 
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3 Existence of self-similar solutions 



In this section we study the steady states of ( 11 ) and ( 12 ) for fixed N and L > N — 6, and pass 



to the hmit e — t- 0. We thus obtain self-similar profiles to the Pradkov model with N < oo that 
turn out to have compact support in [0, A^ — 6]. Afterwards we show that these self-similar 
profiles converge as A^ — > oo. 

3.1 Limit e -^ 



For each e = L/K we choose a steady state g to (11) and (12) with coupling weight T = F"^, 

(24) 



from which we construct a piecewise continuous function gf^ in [0, L] via 

9'n{^k + 0=9n for n = 2...N, k = l...K-l, \C\ < ^e. 



In consistency with the boundary conditions (12) we further define 

gUO = for n = 7...N, < C < ^e, 

as well as 

gf^{L-C) = for n = 2...N, < ^ < |e. 

To ensure that the functions gf^ are well-defined on [0, L] we also set 

gliO = for < ^ < ie, 

but require continuity at ^ = for n = 2...5, that means 

5n(e)=9n(e) for n = 2...5, < ^ < ^e. 



In consistency with (24) we furthermore write 

N 



N 



^^ = ^y„^ p^ = ^(n-6)y,^ 



where 



n=2 



N-6-e/2 



n=2 



N-Q-e/2 



K= I 5^(0 de, Y^= I ^gf^iOd^- 



e/2 



s/2 



This implies 



N-6 



J gUO de = ^n for n = 6...N but J g^O d^ = X', + '-gUO) for n = 2.. .5. 



For the following consideration it is convenient to drop the condition A^ = 1 and to scale the 



steady state solutions differently. Specifically, recalling Remark [9j and due to (23), we can 
assume that 



N 



^{6-n)gf,{0) = ^X^^ = l. 



(25) 



n=2 n=2 

This normalization gives rise to the following uniform BV-estimates. 



12 



Lemma 11. For each < (^ < 1 and M > 6 there exists a constant Cs.m that is independent 
of N, K , and L such that 



' n—6—5 



^=2 n X 71=7 \ n ^ a \ x I 



M, 



n-Q+S 



holds for all N > M and < e < 5, where the measure \d(^g^\ d^ denotes the total variation 
of9n- 



Proof. The assertion follows from the stationary version of (11) since (25) provides uniform 
L^ bounds for g^ and because F"^ is uniformly bounded from above, see Lemma 10 D 



exists a subsequence e — >• such that 



From Lemma 10| Lemma 11| and the normalization condition (|25|) we now infer that there 

(26) 



F^ ^^ r. 



gm 



E^O 



and 



9l 



e-s>0 



£-5-0 



> Qn for n = 2. ..5, 



for n = 2. ..5, 



9n + rUnSn^e for n = 6...N, 



weakly-* in the space of Radon measures M ([0, L]) . Here 6n-6 denotes the delta distribution 
in ^ = n — 6, mQ...mf^ are some nonnegative numbers, and each function gn is nonnegative 



and integrable in [0,L]. It readily follows from (25) and (26) that 

5 



^{6 - n)gn = 1. 



(27) 



n=2 



We now exploit the weak formulation of the stationary version of (11) and show that the 



functions gn satisfy - outside the set of possible singularities - the ordinary differential equa- 
tion ph for self-similar profiles. Moreover, using the weak formulation we also recover the 
boundary conditions and derive algebraic relations for the possible singularities. 

Lemma 12. For each n = 2...N, the function gn satisfies 

L 



j 9n{{i + ^ 



n)(j)' - (/)) d^ + (6 - n)+gn(t){0) - Xn>&nin(t){n - 6) 

= T{j{Jg)^<Pdi + {M4>))^ 

for all smooth test functions (j), where ujn{4>) '■= Xn>6"^n'/'('^ — 6). Moreover, we have 

1. gn G C(/n), 

2. gn has left and right limits at ^ = 0, ^ = n — 7, ^ = n — 5 and ^ = L, 

3. gn G Q}{in), 

4- gn satisfies (JTl) pointwise in In, 



(28) 
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where In '■= (0, L) \ {n — 7 , n — 6, n — 5} and In '■= (0, L) \ {n — 8,n — 7, n — 6,n — 5,n — 4}. 
Proof. We employ Lemma [2] as follows. For a given smooth test function (/> : M — t- M we set 



c , 1 



M^ = /i^ = e-i / md^ = H^k) + 0{e^), 



which gives 






ln = \ f 0(O + (e + 6-n)(A'(Ode + O(e) 



for /c = l...i^-l, 



, 1 



where all error terms depend only on (j) and N. Using (13), and thanks to our definition of 
gniO ^iid gn{L — ^^ ^ C ^ £/2, we therefore find 

L L 

j 9n{{i + Q- nW - 0) de + (6 - n)+ <7n(0)</'(0) = ^ | ( J/)„0de + 0(e). 



The limit e — )• now yields (28). 

Since we have (/>(0) = ijJn{4>) = and J{(jj{(j)))^ = for each (p with compact support in /„, 
we find 

((e + 6-n)5„(e))' = <7n(0+r(J9)„(e) in P'(/„) . (29) 



Since the right-hand side of \2y\ is integrable in /„, we conclude that the function hn defined 
by hn{i) = (C + 6 - n)gn{0 for ^ e (0, L) belongs to W^'^{In) C C(/„), and this implies that 
hn G posesses well-defined one-sided limits at In\ In- Consequently, gn is continuous in /„ 
and has well defined one-sided limits at all points ^ G {0,n — 7, n — 5,L}. Finally, using again 
(29) we deduce that gn is continuously differentiable in /„ and satisfies ph pointwise in that 



set. D 

Lemma 13. The following assertions are satisfied: 

1. We have 

{TKn - l)m„ = for n = 6...iV, (30) 

where Kn denotes the modulus of the n diagonal element of the coupling matrix J, that 
means 

2/3 for n = 2, 

{2p + l)n for n = 3...N-l, 

{I3 + 1)N for n = N. 

2. (a) limgx^o 5n(0) = (jn for n = 2. .A, 
(h) lim^x^o 55(0 = 55 - 6r(/3 -M)m6, 

(c) lim^\^o57(0 = 6r/3m6, 

(d) lim^x^o gnii) =Oforn = 8...N. 
Moreover, lim^ y^ i gn{0 = for all n = 2...N. 
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3. We have 



l[ffn]|(n-7) = r/3(n-l)m„_ 



for n = 8...N 



and 



-|bn]|(n-5) = r(/3 + l)(n + l)m„+i for n = 6...N, 

where \[gn]\{0 ■= lims\o {9n{^ + s) - QnH - «)) ■ 
Proof. 1. Let n = 6...N. Since Qn is integrable in (0, L), we have 

lim(n — 6 lb s)gn{n — 6 it s) =0 
s\0 



(31) 



(32) 



because otherwise Qn would not be integrable. Integrating by parts in (28), and using ([7|, we 
therefore find 

—mn4>{n — 6) = —TKnmn(j){n — 6) 



for all test functions (p that have support in (n — 7, n — 5). Thus we have shown (30). 

2. Let (j) be an arbitrary test function with support in (—1, 1). Combining (28) with luh 
implies, again employing integration by parts and the continuity properties of (?„, that 



lim C5'n(0 + (6 - n)+ gn - S^rUn 1 4>{0) = 
r(6(/3 + l)6l - ^6^6 + 6/3(51) m60(O) 



From this identity we readily derive the claimed formulas for limg\^o ffn (0 ■ Moreover, con- 
sidering test functions (f) with support in (L — 1, L + 1) we find lira^y^i^ §„{£,) = for all 
n = 2...N. 

3. Now let n > 8 and suppose that (j) is supported in (n — 8, n — 6). From (28) and ([7]) 
we now derive 



-\[gn]\{n - 7)</)(n - 7) = r/3(n - l)m„_i0(n - 7), 



which implies (31). The proof of (32) is analogous. 
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As an easy consequence we obtain that all functions gn vanish for ^ > A^ — 6. 

Lemma 14. We have gn{S,) = for all ^ G {N — 6, L] and n = 2...N . 

Proof. Standard results from the theory of ordinary differential equations imply that the 
initial value problem for the system ([7]) with prescribed data at ,^ = L is well-posed on the 
interval (A^ — 6, L\. The claim therefore follows from the boundary conditions for ^ = L, see 
Lemma [Tsl D 



From ( 30 ) we conclude that at most one of the weights run of the Dirac masses does not 



vanish. However, in order to show that all weights vanish we need a better understanding of 
the properties of gn- 
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3.2 Self-similar solutions for N < oo 

In this section we characterise the properties of the functions g„ in more detail, in particular 
the behaviour near ^ = n — 6. The results allow us to conclude that all weights m„ must 
vanish and that the functions gn therefore provide in fact a self-similar solution to the Fradkov 
model. 

All subsequent considerations rely on the solution formula 

gn{n -6±s) = s-2+r«n j t^-^-"g^{n -Q±t)+ f y^-^'^-Gnin - 6 ± y) dy J , (33) 



which is direct consequence of Lemma 12 and the Variation of Constants Principle. Here, 

Gn :=r \B{n- l){gn-i +m„„i5„,_7) + (/? + l)(n + l)(5„+i + m„+i5„_5) j , (34) 
where we set gi = gN+i = and mi = ... = m^ = m^^i = to simply the notation. Notice 



that (33) holds for all n = 2...N and < s < t provided that all terms are well defined, i.e.. 



as long as n — 6 lb s and n — 6 it t belong to [0, L] \ {n — 7, n — 6, n — 5}. 

In order to show that the functions gn are positive almost everywhere on [0, A^ — 6], we 
formulate the following auxiliary result. 

Lemma 15. Suppose there exist 2 < m < N and some ^ £ (0, A^ — 6)\N such that g^iC) = 0- 
Then we have gnU) = for all n = 2...N and all ^ G (0, A^ - 6) \ N. 

Proof. Due to gm > 0, the point ^ is minimizer of gm and thus we have gmiO = OmiO = 0- 
Since gm solves ([T]) pointwise at ^, see Lemma 12, we also find 

5m-i(C) = 9m+l{^ = 0. 

Iterating this argument with respect to n we finally get gn(C) =0 for ^^ i^ = 2...N. In 
particular, we have proven the implication 

9n{C) =0 for some n>2 =^ 9n{^ =0 for all n > 2. (35) 



Since (33) and G2 > imply 

= 52© >(C" + 4) -'+^"^(C + 4)2-^-^92(0 for ^<^<N-6, 



we conclude that 52(C) = for ^ € [^, N — 6). Combining this with (35) we then conclude 
that gniO = for ah ^ G [^, A^ - 6) \ N and n = 2...A^. 

We now choose an index n G N with 7 < n < A^ — 6 such that ^ < n — 6. The solution 



formula ( 33 ) now implies 

= 9n{^>{n-6- C)"'^^"^(n - 6 - e)'"^"'"'5n(0 for < ^ < ^ , 

and arguing as before we derive gn{0 — for all ^ G (0, ^] \ N and n = 2...N. D 

Lemma 16. For each n = 6...N we have Fk^ > 1 and rUn = 0. Consequently, for each 
n = 2...N the function gn is positive and continuous on (0, A^ — 6) \ {n — 6}, and continuously 
differentiable on (0, N — 6) \ {n — 7 , n — 6, n — 5} . 

16 



Proof. Assume for contradiction that there are m > 2 and ^ G (0, A^ — 6) \ N such that 
9m{i) = 0. According to Lemma 15 we have gn{C) = for aU n = 2...N and ^ G (0, A^ — 6) \N. 
It then fohows from Lemma 13 that rriQ = and Qn = for 2 < n < 5, thereby contradicting 



the normahzation condition (27). Consequently, we have g-niO > fo'^ aU n > 2 and ^ € 



(0, A^-6)\N, and the solution formula ([33]) ensures that gniO > for ah ^ G (0, Af-6)\{n-6}. 
Now let n = 7...A^ — 1 be given. Since Gn is nonnegative and gn positive in (0, A^ — 6) \N, 



the solution formula ( 33 ) implies that 

gn{n - 6 ± s) > CnS''^'^^"" for all \s\ < 
where 

mmlgn{n-6- ^),gn{n-6+l)\ 



2' 



i2-rK„ 



>0. 



Since gn is integrable, we now conclude that Tku > 1, and (30) yields m^ 
for rriQ = and rri]^ = are similar. 
Finally, the inclusions 



0. The arguments 



5„ G C((0, A^ - 6) \ {n - 6}) , gn e O{{0, N - 6)\{n - 7,n - 6,n - 5}) 
are implied by = run-i 



m„ 



Corollary 17. We have — < F < 

KQ- - 2-p 



run+i, see Lemma 12 and Lemma 13 
5 



Proof. The upper bound is provided by Lemma [TUl the lower one by Lemma 16 



D 



D 



Since all weights ?tt,„ vanish, we immediately arrive at the following result, which in turn 
implies that the functions gn provide indeed a self-similar profile to the Fradkov model. 

Corollary 18. We have 

1. gn{0) = gn> for all n = 2. ..5 with I]^=2 5n(0) = I, 

2- En=2 Xn = 1 With Xn = J^f "' gn{i) d^ 

3- P = T.n=2in-^)Xn = Q, 

4- A = EI2 Yn > Wtth Yn = /J^~' e5n(6 d^ 

5. T = Fnum/rden depends on (A„)„ and (5'„(0))„ via (|8j). 
Moreover, the weak formulation ^ as well as the identities 

(6-n)5n(0)=X„ + F(JX)„, (6 - n)A„ = F( jy)„ (36) 

hold for all n = 2...N . 

We finally characterize the behaviour of gn near ^ = n — 6. 
Lemma 19. For each n = 6...N one of the following conditions is satisfied: 

1. Fk„ > 2 and gn is continuous at ^ = n — 6 with lim.^^n-6 9n{0 — Gn{n — 6)/(Fk„ — 2), 

2. Fk„ = 2 and gn{n — 6 ± s) ~ —Gn{n — 6) Ins as s — ?■ 0, 
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3. TKn < 2 and gn{n — 6 it s) ~ inS^'^ " for some constant £„ > as s — )• 0. 
Here Gn is defined in ( |34[ ). 
Proof. Throughout this proof we assume that < s < i < 1. We also set 

C„(t):=sup{|G„(0-G'„(n-6)| : ^ ^[n - & - t, n - Q + t\^, 

and notice that Cnif) — )• as t — )• since G„ is continuous at ^ = n — 6. 



Case I: Tku > 2. From (33) we infer that 

9n{n-6±s)=[-) 9n{n-6±t) + ^--^il-[-] 



+ s-2+r«n / yi-r«„ I ^^(^ - 6 ± y) - G„(n - 6) ) dy. 



We therefore find 

Mm sup gn{n — 6 ± s) — Gn(re — 6)/(rK„ — 2) 

and the Umit i — t- provides the desired result. 



< 



i Krt. Z 



Case II: Fk^ = 2. The solution formula (33) gives 



gn{n — 6 lb s) = QniiT- — 6 ± t) + Gn{n — 6)(lnt — Ins) + 



Gnin-6±y)-Gnin-d) 



dy, 



and due to G„(n — 6) > 0, see Lemma 16 



lim sup 

s\0 



we estimate 
5„(n - 6 ± s 



+ 1 



< Cnit). 



Gn{n — 6) Ins 
The claimed asymptotic behaviour now follows by letting i — )• 0. 



Case III: TKn < 2. Formula (33) implies 



1 
s'"^""5n(n - 6 ± s) = Qnin - 6 ± 1) + / y^~^''"Gn{n - 6 ± y) dy, 

s 

and we conclude that the right-hand side of the above equality has a positive limit as s — )• 
0. D 



3.3 Exponential decay and estimates for higher moments 

We next prove that the moments Xn decay exponentially with n where the rate is independent 
of N. This gives rise to tightness estimates that enable us to pass to the limit A^ — t- cxd in 



Section 3.4 Introducing 

(n - l)Xn-l 



^n •" 



nX„ 



r:=— ^>1, $(z):=l + r-- 



we readily derive from (36)i the backward recursion formula 

1 



1 



ZN = T 



Zn = ^{Zn+l) 



for n = 6...N — 1 



(37) 



F/3A^' ' ^^ Tpn 

Notice that <I> is strictly increasing and has exactly two fixed points z = 1 and z = t with 
$'(1) > 1 > *&'(t). We therefore find that z = 1 is unstable, whereas z = r is stable and 
attracts all points z > 1. 
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Lemma 20. We have 

r(l-— )<z„<r for n>N, 

where N is the smallest integer larger than 2(1 + 2/3)^6 = 12(1 + 2f3) . 
Proof. For each n G N with nT > 2 we set y„ := r(l — 2/nr) and find 

^/ N 1 ;^/ N 1 1 + 2/3 2 

^[yn+i) - ^^ — yn> ^[yn) - ^^ — yn 



T(5n ^ - '^ ^ r/3n "^ {iVn Tn-2 

A direct computation reveals that the right hand side is nonnegative for nV > 2(1 + /3), and 
in view of F > 1/ kq we conclude that 

< y„ < Hyn+i) - ^ for n>N. 



For N > N we also have yjsf < zn < t, and using (37) as well as the monotonicity of $ we 



readily verify by induction that Vn ^ Zn < t for all n with N < n < N. D 

Corollary 21. There exist positive constants c and C that are independent of N such that 
CT-'^Xfi < nXn < Cn^/^T~''Xfi for N <n<N . 



Proof. By Lemma 20 we have 

m=N+l m=N+l 

The concavity of the logarithm implies In (l — ^) > — Ym-2 ' ^^'^ hence 
, / -A- / 2 \\ /•" ds 2, /Fn-2 

m=Af+l -"^ 

We therefore find 

iVX^7 ^ _.7/FiV-2\2/r 



- nX„ - V Fn - 2 / 



.n-N ^ 



and this implies the desired result since F is bounded, see Corollary 17 D 



We now exploit the exponential decay of X„ and derive tightness estimates. To this end 
we consider the moments 



Jo 



and notice that X„ = Mo^n and Yn = Mi.„. 

Lemma 22. For any k >0 there exists a constant Ck > independent of N such that 

N N 



^n''Xn + ^Mk,n<Ck. 



n=2 n=2 
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Proof. Let k > 1. Multiplying ([7]) by ^ , integrating over (0, A^ — 6) and using integration by 
parts as well as the boundary conditions, we find 

{k - l)Mfc,„ = k{n - 6)Mfc_i,„ + r{JMk)n . 

Summing over n and using (|15|)i we deduce 



N N N 

ik-l)J2 Mk,n = k ^(n - 6)Mfc_i,„ < A; ^ nMfc_i , 

n=2 n=2 n=2 

Since Holder's inequality for integrals implies 



we can employ Holder's inequality for series to find 

N N , 



^ Af (fc-l)/fc. " ^ .1/k 



n=2 n=1 n=2 

and hence 

k \k. 



n=2 n=2 



Thanks to Corollaries 17 18, and 21 we then obtain 



N N 

n.=2 n=N+l 

and this completes the proof for A; > 1. The case k G (0, 1] follows by interpolation. D 

We finally prove that even some moments with exponential weight are uniformly bounded. 

Lemma 23. For each < A < Inr there exists a constant C\ that is independent of N such 
that 

N N ^ 

n=2 n=2"^'^ 

Proof We multiply ([T]) by e"^^ and integrate over (0, A^ — 6) to obtain 

rN-6 

(6-n)5„(0)+ / e^«(A(e + 6-n)-l)5„(e)de = r(Ji?A),, 

Jo 

where -Ea,^ = /q e^^g'„(^) d^, and this implies 

rN~6 rN-6 

Jo Jo 



Now we choose A with < A < A and A^ < A Inr, and estimate 

-Af-6 _ ^ f\n/X ^ ^ rN-G 

lxn/\ 



/ ne^^gniOdC= ne""^ QniO dC + ne^^gniOdC 

Jo Jo JXn/\ 
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This implies 



N~6 



(A - A) / (e^^gniO d^ < Ane^ "/^X„ + r{JEx)^ 
Jo 



and summation over n > 2 yields, thanks to (15)i, 



N 



(A - A) j; / ^e'^gniO ^^ < A J] e''^/'X„ 

n=2"'0 n=2 



Finally, the estimate 



N 



N 



5^e^"X„ + J]e^VAj^„<C7A 



n=2 



n=2 



TV 



follows from Corollary 21 due to X^ < X^,„=2 "^^n^ ~ ^ ^^^ ^^^ choice of A and A. 
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3.4 Limit A^ -^ oo 

To finish our existence proof we construct self-similar profiles to the original Fradkov model 
by passing to the limit N ^ oo. Our arguments are very similar to those used in Sections 3.1 

0, and therefore we only sketch the main ideas. 

[0, oo) 



and 3.2 for the limit e 



|0, oo) with n > 2 as trivial 



For each A^ < oo we define the functions g^ 
continuation of the self-similar profiles constructed in Section 3.2 This means we set g^iO — 
for n>iVor^>A^ — 6. We also denote the corresponding coupling weights by P and use 
notations such as X^ ,YJ^ and M^^ to refer to the various moments of g^ . 



Due to the BV-estimates from Lemma 11 and the normalization, see Corollary |18[ there ex- 
ist a subsequence N — )■ oo, nonnegative real numbers {j)n)2<n<b and ("T'n)„>e5 and a sequence 
of nonnegative integrable functions {gn)n>2 such that 



ffn(O) > 



for n = 2. ..5, 



and 



/V N—^oo 
9n ^ 



gn + Xn>6i^nSn-6 weakly-:*r in M([0, do)) for n > 2. 



Moreover, the uniform moment estimates from Lemma 22 imply gn G L^((0, oo);^'^(ii^) for all 
A; > 0, and hence 



z2^^ 



N- 



^Xn. 



E^^n 



7V->oo 

k 



y nXn < oo 



n=2 



with Xn = rUn + f^ gniO d^- We now conclude that F^ — )• F, where F depends self- 
consistently on the boundary data gn for 2 < n < 5, and the moments Xn- 



We are now in the same situation as in Sections 3.1 and |3.2[ In particular, analogously 
to the proofs of Lemma 12 and Lemma 13 we show gn = 5n(0) for 2 < n < 5, mn = for 



n > 6. We then infer that {gn)n>2 provides a weak self-similar solution to the Fradkov model 
with N = oo, which satisfies the assertions of Corollary 18 and Lemma 19, Finally, it is clear 



by construction that this self-similar solution satisfies the moment estimates from Lemma 23 
and thus we have finished the proof of Theorem [T| 
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4 Numerical examples 



To illustrate our analytical results we implemented the explicit Euler scheme for (11), which 



has the notable property that all moment balances, see Lemma [2] and Corollary [4j remain 
valid provided that we replace the continuous time derivative by its discrete counterpart. 
In particular, computing F by ( |18[ ) our scheme conserves the area and polyhedral defect 
up to computational accuracy. Moreover, due to the upwind discretization of the transport 
operators, and since Lemma [6] provides upper and lower bounds for F, one easily shows that 
the explicit Euler scheme preserves the nonnegativity of the data provided that the time step 
size is sufficiently small. 
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Figure 4: Self-similar solution for several values of /3, where circles and lines represent numer- 
ical data and interpolating splines, respectively. To-p row: F and "^^=2-^^ versus /?. Bottom 
row: Xn and Yn versus n with vertical line at n = 6. 



We performed a large number of numerical simulations for various values of /3 and different 
types of initial values (random, uniformly distributed, and several variants of localized data). 
To ensure that the initial data belong in fact to the set Z//_|_, we chose at first nonnegative 
values g^ for n = 2...N — 1 and k = 1...K — 1, and computed afterwards two scaling factors 
Qi and 02 such that g^ = (x2<n<5ai + X6<n<NCt2)gn satisfy the constraint P + eQ = and 
yield the prescribed area. 

In our simulations we observed that all numerical solutions for a given value of (3 converge, 
as t — )• oo, to the same steady state. We therefore conjecture, that for all (3 and N < oo there 



exists a unique steady state that is moreover a global attr actor for (11). It would be highly 



desirable to give a rigorous justification for this numerical observation, but even to prove the 
uniqueness of F remains a challenging task. We also conjecture that for N = oo there is only 
one self-similar solution with fast decay in ^, but emphasize that self-similar solutions with slow 
decay might exist as well. Such solutions exist in related mean-field models for coarsening that 
couple transport and coalescence [13], but cannot be detected by our approximation scheme. 

The numerically computed steady states for several values of /3 € [0.2, 5.0] are, along with 
some derived data, depicted in Figures [ij [Sj and [6) All computations are performed with 
A = 1, N = 25, L = 20, and e = 0.05, and due to the numerically computed residuals we 
expect that the discrete solutions g^ resemble the limit profiles gn with N = oo very well. In 



22 



Qn versus g for n=2 



+3.1E-2 




Qn versus ^ for n=4 





3. 1 


2. 3. 4. 

Qn versus f for n=6 


+ 1.0E0 
5.0E-1 


V 


^V^ 


0. 





+3.1E-3 




+ 1.2E-4 



2. 4. 6. 

ffn versus f for n=20 




+ 1.2E-1 



+5.4E-2 



+9.1E-4 




+4.6E-5 



15. 20. 



5. 10. 15. 20. 

Qn versus g for n=25 




15. 20. 



Figure 5: Plots of Qn versus ^ with vertical lines at .^ = n — 5, .^ = n — 6, and ^ = n — 7. Notice 
that the plot range for ^ varies with n. 
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Figure 6: Decay of the solutions: In ( Yln=2 5"(0) versus ^ and lnX„ versus 



n. 



particular, Figure [6] confirms that the self-similar profiles for N = oo decay exponentially in 

Figure Isl illustrates that for n > 7 there is no pointwise convergence g^ > gn at the 

critical point ^ = n — 6. In fact, at least for small /3 and moderate values of n we observe 
that the discrete data g^" are considerably smaller than gn{n — 6). This phenomenon stems 
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from our discretization and can be understood as follows. On the discrete level steady states 



satisfy, see (11) 



{TKn - l)5^ = r(/3(n - 1)5^1 + (/3 + l)(n + l)gll^ , 
{TKn - l)gt^' - gt^' = r(/3(n - l)gtf + (/? + l)(n + l)gtlf 



For small e we can express the right hand sides in terms of gn±i (n — 6), the values of the limit 
functions gn±i at ^ = n — 6. Equating the resulting right hand sides we then conclude that 

{TKn - l)gt = (r^n - l)9t^' - gt^' + oil). 



If the limit function gn is also continuous at ^ = n — 6 (according to Lemma 19, this happens 
for r > 2/Kn = 2/(2/3 + l)n and hence at least for large n) we can approximate the terms 
gn"^^ and g^"^^ by ^^(n — 6). This gives 

""" ^ ^.~(~-e) = (i- ,p^;,,„_, >„(n-o). 

where the right hand side is always nonnegative due to 5n(n — 6) > and Tku > 2. 
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